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The  dynamics  of  single-crystal  silicon  —100  jam  size  rectangular  paddle  oscillators  at  room 
temperature  have  been  studied  using  a  recently  developed  high-resolution  scanning  laser  vibrometer. 

The  dynamic  mechanical  behavior  is  determined  by  scans  of  the  entire  device,  providing  both 
amplitude  and  phase  spatial  maps  of  the  vibratory  response.  These  reveal  more  than  16  normal 
modes  below  500  kHz.  In  addition  to  simple  translation  and  torsional  motion,  flexural  modes  of  the 
paddle  plate  are  observed.  Quality  factors  ranging  from  1  X  103  to  2X  104  are  measured  and  are 
found  to  be  significantly  lower  than  those  expected  from  well-known  intrinsic  absorption 
mechanisms.  The  measurements  reveal  that  there  exists  significant  modification  of  the  expected 
eigenfrequencies  and  mode  shapes.  It  is  speculated  that  this  is  caused  by  excessive  undercutting  of 
the  support  structure,  and  that  the  resulting  energy  flow  into  the  support  leads  to  increased  oscillator 
loss.  Indeed,  some  correlation  is  found  between  observed  loss  and  energy  levels  resident  in  the 
supports.  At  frequencies  where  there  is  relatively  little  support  motion,  three-dimensional 
finite-element  modeling  accurately  predicts  the  paddle  modal  behavior.  ©  2002  American  Institute 
of  Physics.  [DOI:  10.1063/1.1502014] 


I.  INTRODUCTION 

High- Q  mechanical  oscillators  based  on  single-crystal 
silicon  microfabrication  technology  are  of  great  interest  in  a 
variety  of  applications.  Examples  include  charge  detectors,1 
radio-frequency  filters,2  magnetic  resonance  force 
microscopes,3  and  torque  magnetometers.4  Up  to  the  present 
time,  ultra-high-2  (106-108)  have  been  observed  only  in 
oscillators  having  at  least  centimeter  lateral  dimensions,  a 
good  example  of  which  is  the  so-called  “double  paddle” 
type.5-7  In  microelectromechanical  systems  (MEMS),  how¬ 
ever,  2  are  found  to  be  lower  than  expected  based  on  scaling 
considerations  of  fundamental  loss  mechanisms.8  Indeed,  re¬ 
cent  measurements  by  Evoy  et  al.9  made  on  5. 5 -/zm- square 
Si  paddles  demonstrated  Q  reaching  only  2.5  X  103  at  room 
temperature. 

The  specific  cause  of  low  Q  in  MEMS  oscillators  is  not 
presently  understood.  At  larger  length  scales,  experiments10 
have  shown  that  eliminating  metal  films  used  for  capacitive 
driving  and  detection  may  increase  Q  as  much  as  threefold. 
Recent  experimental  results  on  macroscale  single-crystal  sili¬ 
con  double-paddle  oscillators  suggest  that  in  these  high- 2 
oscillators  the  internal  friction  in  lower- 2  modes  can  be  cor¬ 
related  with  vibration  of  the  supporting  structure.11  Indeed,  it 
is  for  those  modes  with  minimal  coupling  to  the  supports  that 
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moderately  high  Q  have  been  obtained  (Q~  106  at  300  K).  In 
nanofabricated  single-crystal  silicon  wires,  Q  has  been 
shown  to  decrease  as  the  surf ace-to- volume  ratio  increases, 
suggesting  that  surface  loss  mechanisms  play  an  increasing 
role  as  oscillator  size  is  reduced  to  the  nanoscale.12  To  im¬ 
prove  our  understanding  of  the  loss  mechanisms  at  work  in 
silicon-based  micro-oscillators  in  order  to  achieve  higher- 2 
versions  of  MEMS  oscillators  in  general,  we  have  investi¬ 
gated  the  dynamics  and  loss  mechanisms  of  some  MEMS 
oscillators  with  a  recently  developed,  high-resolution  in  situ 
measurement  technique  capable  of  revealing  the  detailed  mo¬ 
tion  of  both  the  oscillator  and  the  supporting  structures  as 
well.  In  this  article,  we  report  such  high-resolution  (<3  jam) 
laser  Doppler  vibrometry  scanning  measurements  made  on 
single-crystal  silicon  paddle  oscillators  having  lateral  spatial 
dimensions  in  the  100  jam  range.  The  resulting  detailed  in¬ 
formation  is  used  to  identify  the  actual  modal  response  of  the 
complete  structure  in  an  attempt  to  understand  the  low- 2 
values  seen  in  general  for  such  microscillators. 

II.  SAMPLE  PREPARATION 

The  oscillators  were  fabricated  on  a  silicon-on-insulator 
(SOI)  (100)  substrate  using  optical  lithography  and  Si  micro¬ 
machining  techniques  at  the  Cornell  Nanofabrication 
Facility.13  The  top  Si  layer,  which  after  processing  makes  up 
the  oscillator  and  the  top  layer  of  the  supporting  structures,  is 
1.5  jam  thick.  The  oxide  layer,  which  is  removed  to  produce 
released  features,  is  3  jam  thick.  The  silicon  is  nominally 
undoped  with  a  resistivity  of  14-22  and  1-3  12  cm  for  the 
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FIG.  1.  Scanning  electron  micrograph  of  a  single-crystal  silicon  paddle  os¬ 
cillator  fabricated  in  SOI  substrate  with  d=  147  jmm.  An  array  of  eight  os¬ 
cillators  with  d=  100-270  jmm  was  fabricated  on  a  single  wafer.  Other  pa¬ 
rameters  are  (in  /un):  w  =  97.8,  L  =  50.9,  b  =  4.15,  S{=  130.0,  S2  =  183.8, 
and  S3  =  23.5;  all  with  thickness  h=  1.5.  The  oxide  layer  etched  to  release 
the  oscillator  was  3  jmm. 

top  Si  and  substrate,  respectively.  Before  processing,  the  wa¬ 
fer  was  rotated  45°  with  respect  to  the  flat  so  that  the  arms  of 
the  oscillator  are  parallel  to  the  100  direction.  Eight  such 
oscillators  of  differing  lengths  were  fabricated  on  the  same 
wafer.  For  the  eight  oscillators,  the  paddle  length,  d,  varies 
from  97.6  to  268  jmm  (Fig.  1)  while  all  other  dimensions 
remain  fixed:  (in  jmm)  w  =  97.8,  L  =  50.9,  b  =  4.15,  sq 
=  130.3,  s2=  183.8,  and  ^3  =  23.5.  This  study  was  focused  on 
the  oscillator  having  d=147  jmm.  Figure  1  is  a  scanning 
electron  micrograph  of  this  oscillator.  The  rectangular  paddle 
plate  was  designed  with  5-yam-square  perforations  to  assist 
in  the  etching  of  the  oxide  layer  under  the  large  area  of  the 
plate. 

III.  MEASUREMENT  APPARATUS 

The  scanning  laser  Doppler  vibrometer  microscope 
(FDVM)  instrument  we  developed  for  this  study  is  com¬ 
prised  of  three  main  subsystems:  The  optical  interferometer, 
the  demodulator  and  recording  electronics,  and  a  three-axis 
mechanical  scanner.  A  schematic  of  the  system  is  shown  in 
Fig.  2.  The  optical  interferometer  system  is  a  quadrature 
heterodyne14  laser  interferometer  assembled  from  both  bulk 
and  fiber-optic  components.  It  incorporates  an  argon-ion  la¬ 
ser  (Spectra  Physics  model  2060-4  with  a  514.5  nm  etalon 
for  single-frequency  operation).  A  variable  half-wave  plate 
(Fambda  Research  Optics)  is  employed  at  the  laser  output  to 
control  the  balance  of  power  in  each  arm  of  the  interferom¬ 
eter.  The  beam  then  passes  through  a  polarizing  beam  splitter 
(PBS)  (Fambda  Research  Optics,  PB-12B-R-488-514).  The 
reflected  component  from  the  PBS  is  vertically  polarized  and 
launched  down  a  polarization  maintaining  fiber  (Alcoa 
Fujikura  SM48-P-S,  450  nm  cutoff).  The  fiber  is  used  to 
deliver  light  to  a  compact  probe  housing  containing  a  pair  of 
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FIG.  2.  Schematic  diagram  of  the  scanning  laser  Doppler  vibrometer  micro¬ 
scope  (LDVM)  system  and  associated  demodulation  setup.  Notation:  PBS, 
polarizing  beam  splitter;  PM,  polarization  maintaining  fiber;  PD,  photode¬ 
tector;  /X  2,  frequency  doubler;  BPF,  band  pass  filter;  SAW,  surface  acous¬ 
tic  wave;  and  LPF,  low  pass  filter. 

axial  GRIN  lenses  (Newport  Optics,  GPX-043  AR  14)  sepa¬ 
rated  by  a  quarter-wave  plate  (Fambda  Research  Optics,  WP- 
12-CQ-M-4-514.5).  The  light  then  passes  through  the  optical 
window  of  a  vacuum  chamber  to  the  sample  surface.  For  the 
measurements  reported  here,  the  typical  spot  size  is  2.5  jmm 
in  diameter,  and  the  typical  power  levels  are  kept  below  0. 1 
mW  in  order  to  insure  minimal  sample  heating.  On  the  return 
path,  the  reflected  light  from  the  sample  passes  through  the 
GRIN  lenses  and  quarter-wave  plate  assembly,  thereby 
achieving  a  rotation  of  90°  relative  to  the  input.  This  rotation 
insures  horizontal  polarization  and  maximum  transmission 
through  the  PBS  to  the  photodetector.  The  GRIN  lenses 
minimize  spherical  aberration  so  that  a  diffraction-limited 
spot  can  be  achieved.  These  optical  elements  along  with  the 
optical  fiber  are  assembled  into  a  compact  housing  mounted 
to  a  three-axis  positioning  stage. 

The  transmitted  beam  through  the  PBS  comprises  the 
reference  beam  of  the  interferometer.  It  passes  through  a  35 
MHz  Bragg  cell  (Isomet,  Acouto-Optic  modulator  model 
120 IE)  and  quarter-wave  plate  (Fambda  Research  Optics, 
WP-12-CQ-M-4-514.5),  where  it  is  rotated  45°  relative  to 
the  horizontally  polarized  incident  light.  This  is  done  in  order 
to  impose  circular  polarization  on  the  reference  beam.  A  mir¬ 
ror  is  aligned  to  retroreflect  the  first-order  diffracted  beam 
back  along  this  path  to  the  Bragg  cell  and  the  PBS.  As  a 
result,  the  polarization  rotation  and  frequency  modulation  of 
the  reference  beam  take  place  twice;  thus,  the  reference  beam 
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is  up-shifted  70  MHz  and  rotated  90°  relative  to  the  incident 
beam. 

Both  the  frequency-modulated  reference  beam  and  the 
reflected  signal  from  the  sample  surface  pass  back  through 
the  PBS,  a  Gian  Thompson  polarizer  (Melles  Griot 
03PTH001/A),  and  are  collected  at  a  ~1  GHz  avalanche 
photodiode  receiver  (New  Focus  1601  AC).  The  polarizer 
serves  to  insure  a  common  polarization  upon  mixing  of  the 
two  beams. 

A  hybrid  analog -digital  approach  is  employed  for  the 
demodulator  and  recording  subsystem.  The  photodetector 
contains  a  frequency-modulated  (FM)  electrical  signal  carry¬ 
ing  a  phase  signal  that  is  the  dynamic  surface  displacement 
of  the  target  normalized  to  the  wavelength  of  light.  The  pho¬ 
todiode  voltage  can  be  written  as14 

I  d(t)  \ 

ysig=  V  sinl  ( obt  +  +  (4>{t))\  +  Vnoise , 

where  V  is  a  constant,  cob  is  the  frequency  shift  due  to  two 
passes  through  the  Bragg  cell,  d{t)i 47tX  is  the  phase  modu¬ 
lation  due  to  the  out-of-plane  displacement  of  the  sample 
surface  referenced  to  the  wavelength  of  the  laser  (X),  (<f>(t)) 
is  the  laser  phase  noise,  and  Vnoise  is  the  noise  associated 
with  the  photodiode  and  receiver  system.  The  dominant 
source  of  noise  in  this  system  is  associated  with  the  phase 
noise  of  the  laser  and  determines  the  minimum  detectable 
level  (4X10-5  nm/^Hz).  The  maximum  displacement  for 
these  measurements  is  approximately  5  nm.  It  should  be 
pointed  out  that  this  system  is  capable  of  measurements  up  to 
—500  MHz  for  the  sub  wavelength  surface  displacements 
typically  associated  with  MEMS  resonator  devices. 

The  phase-modulated  signal  plus  noise  is  extracted  using 
a  hybrid  analog -digital  method.  The  front  end  of  the  de¬ 
modulator  (Fig.  2)  employs  an  analog  stage  (various  Mini- 
Circuits,  Inc.,  components)  to  extract  the  quadrature  compo¬ 
nents  at  the  base  band  of  the  system.14  These  are  digitized, 
and  simple  processing  is  performed  to  determine  the  arctan¬ 
gent  (i.e.,  phase)  of  the  signal.  This  is  accomplished  by  gen¬ 
erating  a  pair  of  electrical  reference  signals  phase  locked  to 
the  Bragg  cell  drive  signal  and  phase  shifted  90°.  This,  in 
turn,  is  done  by  tapping  off  the  35  MHz  Bragg  cell  electrical 
drive  signal,  doubling  it,  filtering  it,  splitting  it,  and  retarding 
one  of  the  outputs  by  90°  (Fig.  2).  These  two  reference  sig¬ 
nals  are  then  individually  mixed  with  the  output  of  a  splitter 
fed  by  the  FM  signal  from  the  photodetector.  The  two  signals 
are  subsequently  captured  with  a  14  bit  wavelength  digitizer 
(Datel  model  PCI-4 16N).  The  displacement  plus  noise  nor¬ 
malized  to  the  wavelength  of  the  laser  light  is  then  found  as 
the  output  of  a  four-quadrant  arctangent  routine  on  a  com¬ 
puter,  viz.  tan _1(//2)X/4 77^2,  where  I  and  Q  are  the  two 
quadrature  signals  (Fig.  2),  X  is  514.5  nm,  and  n  is  the  index 
of  the  medium  where  the  Doppler  shift  occurred  (n  =  1). 


IV.  EXPERIMENTS 

Scans  were  performed  on  each  of  the  oscillators  at  room 
temperature.  A  mechanical  scanner  subsystem  was  used  to 
precisely  position  the  probe  relative  to  the  sample.  For  these 


experiments,  the  silicon  wafer  sample  was  mounted  with  car¬ 
bon  tape  on  a  500  kHz  ultrasound  transducer  (K-B  Aerotec 
0.5  MHz- Alpha)  that  was  in  turn  mounted  in  the  vacuum 
chamber.  The  vacuum  was  maintained  to  <0.1  mTorr  to 
eliminate  acoustic  radiation  damping.  Scanning  motion  takes 
place  in  air  outside  the  vacuum  chamber  with  the  probe  head 
portion  of  the  vibrometer  mounted  on  a  three-axis  scanning 
stage  (Dadel  8122P-0414)  with  incremental  optical  linear  en¬ 
coder  feedback  with  100  nm  resolution  (Dynamic  Research 
Corporation,  model  LB). 

In  order  to  study  the  motion  of  the  entire  structure,  the 
scan  area  was  made  large  enough  to  include  the  support  re¬ 
gions.  For  the  measurements  reported  here,  the  probe  was 
scanned  over  the  paddle  and  boundary  in  10  jmm  steps.  A  step 
size  of  2  jmm  was  used  for  the  supporting  arms.  The  sample 
locations  on  the  oscillator  plate  region  were  selected  to  miss 
the  perforations.  The  wafer  was  attached  to  a  pieozoelectric 
transducer  crystal  driven  by  a  linearly  swept  frequency  chirp 
from  1  to  500  kHz  with  a  temporal  length  of  4  /xs.  A  typical 
measurement  included  acquisition  of  the  associated  transient 
response  at  3141  locations  with  100  averages  per  point. 
These  records  were  then  windowed  (Hanning  window)  and 
transformed  by  fast  Fourier  transform.  With  respect  to  the 
oscillator  upon  which  we  concentrated,  the  magnitude  and 
phase  maps  for  eight  of  the  more  prominent  resonances  be¬ 
low  500  kHz  are  plotted  in  Fig.  3.  Nodal  lines  appear  as  deep 
blue  and  large  amplitude  regions  as  red. 

Four  of  the  six  lowest- frequency  modes  observed  [Figs. 
3(a),  3(b),  3(c),  and  3(e)]  are  principally  “rigid-body-”  type 
modes  of  the  paddle  plate,  i.e.,  the  plate  remains  planar  while 
the  arms  and  supports  undergo  bending.  Modes  (d)  and  (f)- 
(h)  display  considerable  flexural  plate  motion.  In  particular, 
(f),  (g),  and  (h)  resemble  the  classic  modes  of  a  plate  with 
four  free  boundaries. 


V.  DISCUSSION 

In  addition  to  motion  of  the  oscillator  plate  itself,  the 
scan  readily  reveals  motion  in  the  support  structure.  The  con¬ 
siderable  motion  of  the  top  Si  support  layer  (the  arm  appear¬ 
ing  on  the  left  in  Fig.  3)  suggests  that  the  SOI  underlying 
these  structures  has  been  undercut  by  the  etching  process. 
Scans  of  the  other  oscillators  reveal  similar  mode  structure  as 
in  Fig.  3,  but  with  different  resonance  frequencies,  which  are 
consistent  with  predictions  based  on  the  overall  change  in  the 
plate  dimensions. 

A  finite-element  model,  created  using  the  program 
Sara3D,15  was  employed  to  predict  the  dynamic  behavior 
expected  for  the  oscillator  assuming  it  was  tied  to  a  rigid 
support  structure.  The  model  used  orthotropic  silicon  mate¬ 
rial  properties.16  It  required  4040  solid  elements  and  discreti¬ 
zation  of  approximately  two  elements  per  perforation  and 
two  elements  through  the  thickness.  The  excitation  was  a 
harmonic  displacement  applied  to  the  oscillator  at  one  of  the 
plate  corners.  For  modes  where  the  plate  vibrations  have 
nodes  at  the  attachments  minimizing  coupling  to  the  sup¬ 
ports,  the  calculation  accurately  predicts  both  the  frequency 
and  mode  shape  experimentally  observed.  For  example,  for 
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(e)  (f)  (g)  00 

FIG.  3.  (Color)  LDV  scan  (amplitude,  upper;  phase,  lower)  of  the  paddle  oscillator  and  supporting  structure  decomposed  by  frequency.  The  displacement  of 
8  of  the  19  modes  observed  below  500  kHz  are  shown  with  the  amplitude  color  scale  in  meters  and  the  phase  scale  going  from  red  (+  tt)  to  green  (—  7r).  Modes 
(a),  (b),  (c),  and  (e)  are  rigid-plate-type  modes,  while  (d)  and  (f)-(h)  involve  considerable  plate  bending.  Mode  (d)  is  similar  to  the  fundamental  mode  of  a 
beam  clamped  at  both  ends.  Mode  (h)  exhibits  the  largest  Q  of  all  modes  observed. 


the  [2,2]  mode  of  Fig.  3(h),  the  measured  and  predicted  fre¬ 
quencies  are  within  1%.  For  the  resonances  in  which  the 
oscillator  is  expected  to  have  significant  motion  near  the  sup¬ 
port,  the  modal  behavior  and  frequency  actually  observed 
were  not  predicted  well.  It  is  reasonable  to  assume  that  the 
change  in  mode  structure  and  eigenfrequencies  is,  in  fact, 
due  to  excessive  undercutting  of  the  support  structure. 

Table  I  lists  Q  ~ 1  determined  at  the  frequencies  displayed 
in  Fig.  3  together  with  some  additional  lower-amplitude 
resonances.  Here,  Q  range  from  about  1  X  103  to  2X104. 
Since  our  technique  provides  detailed  spatial  maps  of  the 
vibratory  response,  we  are  able  to  determine  other  useful 
oscillator  parameters.  For  example,  also  shown  in  Table  I  is 
the  kinetic  energy  in  the  bottom  support  and  the  attachment 
structure  normalized  to  the  total  kinetic  energy  as  determined 
from  the  LDVM  scans  by  integrating  the  square  of  the  ve¬ 
locity  ( v 2)  over  the  pertinent  structures.  This  metric  allows 
us  to  test  for  a  quantitative  correlation  between  the  amount 
of  kinetic  energy  near  the  attachment  points  to  the  base  struc¬ 
ture  and  the  system  loss11  ( Q~l ).  One  observation  that  can 


be  made  in  Table  I  is  that  for  all  the  modes  having  Q 
>  104,  the  energy  ratio  is  quite  low,  i.e.,  no  more  than  several 
tenths  of  1%.  Conversely,  the  very  lowest  Q  (those  ~103) 
have  relatively  high  ratios,  i.e.,  on  the  order  of  10%.  In  be¬ 
tween  these  high  and  low  values  of  Q ,  there  is  a  less  clear 
correlation  of  Q  and  the  energy  ratio. 

The  known  intrinsic  absorption  mechanisms  for  crystal¬ 
line  silicon  do  not  account  for  these  measured  values  of  loss 
even  for  those  modes  exhibiting  the  highest  Q.  Specifically, 
for  phonon-phonon  interaction17  caused  by  crystal  lattice 
anharmonicity  one  would  expect  at  frequency  /  to  observe 
/2~1014  Hz;  for  thermoelastic  dissipation  induced  by  lon¬ 
gitudinal  inhomogeneous  dynamic  strains18  one  would  an¬ 
ticipate  /2~1016  Hz;  and  for  thermoelastic  dissipation  as¬ 
sociated  with  transverse  strains19  fQ~  1012  Hz  is  expected. 
We,  therefore,  predict  that  these  internal  dissipation  mecha¬ 
nisms  would  result  in  much  lower  losses  than  those  shown  in 
Table  I.  Indeed,  our  highest  fQ  product  is  found  for  the  high¬ 
est  frequency,  lowest  loss  mode  and  is  fQ  =  8.5X109  Hz. 
This  is  two  orders  of  magnitude  below  the  limit  set  by  ther- 
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TABLE  I.  Mode  frequencies,  measured  Q' s,  and  energy  ratio  ( v 2  level  over 
the  support  structure  divided  by  v2  over  entire  structure). 


Mode  frequencies 

0X1O“3 

Energy  ratio  X  102 

33  600 

10.6 

0.00099 

48  645 

17.3 

0.00061 

62.836 

1.1 

11.6 

72  021 

1.2 

11.3 

82  916 

1.2 

12.7 

87  891 

2.7 

4.9 

115  021 

14.2 

0.0015 

150  665 

2.7 

21 

189  056 

2.1 

17.2 

222  656 

1.7 

0.72 

310516 

12.5 

0.17 

373  383 

9.0 

0.016 

419  647 

11 

0.061 

441  315 

9.0 

1.2 

456  360 

18.6 

0.0033 

moelastic  dissipation.  We  thus  speculate  that  for  all  the 
modes  of  this  oscillator,  the  observed  Q  is  in  fact  controlled 
by  energy  flow  through  the  undercut  supports  and  through 
the  attachments,  with  the  highest-<2  modes  being  those  that 
are  most  isolated  from  the  supports  because  of  their  mode 
shape.  For  the  data  reported  here  the  [2,2]  mode  at  456  kHz 
has  the  lowest  measured  loss,  which  is  consistent  with  the 
small  force  moment  expected  for  this  mode  at  the  attachment 
of  the  plate  to  the  support  [Fig.  3(h)].  It  is  anticipated  that  a 
modification  of  the  design  of  these  devices  that  both  immo¬ 
bilizes  the  support  structure  by  eliminating  undercutting  and 
also  minimizes  the  transfer  of  energy  to  the  attachments 
could  result  in  the  observation  of  Q  controlled  by  intrinsic 
thermal  elastic  dissipation.  This  should  result  in  an  increase 
in  the  measured  Q  by  as  much  as  two  orders  of  magnitude. 
Work  is  in  progress  focused  on  these  goals. 
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